Optical pulse propagation in a switched-on photonic lattice: Rabi effect with the roles 

of light and matter interchanged 
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A light pulse propagating in a suddenly switched on photonic lattice, when the central frequency 
lies in the photonic band gap, is an analog of the Rabi model where the two-level system is the 
two resonant (i.e. Bragg-coupled) Fourier modes of the pulse, while the photonic lattice serves as 
a monochromatic external field. A simple theory of these Rabi oscillations is given and confirmed 
by the numerical solution of the corresponding Maxwell equations. This is a direct, i.e. temporal, 
analog of the Rabi effect, additionally to the spatial analog in the optical beam propagation described 
in Opt. Lett. 32, 1920 (2007). An additional high-frequency modulation of the Rabi oscillations 
refiects the lattice-induced energy transfer between the electric and magnetic fields of the pulse. 



The propagation of classical waves in periodic struc- 
tures has been known for a long time to exhibit intrigu- 
ing analogies of the quantum phenomena, such as Bloch 
oscillations [ij and Zener tunneling [2]. Optical demon- 
stration of these two effects have been performed in the 
one-dimensional periodic structures based on waveguide 
arrays and superlattices d, 0J1, 0, S B Q , and recently 
in the two-dimensional case [10| . In this letter it is shown 
that there is yet another type of analogy in the propa- 
gation of optical pulses in the periodic photonic lattices, 
where the roles played by the light and matter are inter- 
changed. 

An optical pulse with the central frequency lying inside 
the photonic band gap is Bragg reflected by the lattice, 
which changes its spatial Fourier index by the recipro- 
cal lattice vector. Applying the same argument to the 
reflected pulse, we see that there is a resonant coupling 
between the two optical modes with the spatial Fourier 
indices lying on the opposite sides of the Brillouin zone. 
The amplitudes of the resonant Fourier modes are sub- 
ject to oscillations induced by the periodic lattice, which 
can be called Rabi oscillations where the optical pulse 
(its two resonant modes) plays the role of a two-level 
system and the photonic lattice takes the place of an ex- 
ternal monochromatic field. Below we provide a simple 
theory of this effect and confirm the predictions by direct 
numerical simulations. 

The spatial analogs of Rabi oscillations in an optical 
beam propaga ting in a photonic crystal have been previ- 
ously studied [Til [l^ (see also Ref. The temporal 
oscillations studied here is a direct analog of Rabi effect. 

The optical pulse propagation (oblique, in general, see 
Fig. [T|) is described by the following equations for the 
TE-like and TH-like pulses [H: 
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where E = e^E^^^') and H = e^i?(™), and V 
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the TM-like case can be treated similarly. 

The refraction index s{x) consists of the uniform com- 
ponent So and a periodic modulation (a weak photonic 
lattice), i.e. 
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where kB = i^/d with d being the lattice period and 
the Fourier amplitudes of the lattice satisfy «_„ = w* 
(below, the uniform refraction index Eq is accounted for 
by introducing a modified speed of light cq = c/ 
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eydy. Below we will concentrate on the TE-like pulses. 
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FIG. 1: (Color online) Schematic setup. The electric and 
magnetic fields for the TE-like and TM-like pulse propagation 
are shown by the arrows, the photonic lattice is represented 
by the bars. 

Equation ([T]) for the TE-like pulses can be considered 
in Fourier space by setting E{x,y,t) — ^'(a;, t)e*''^ and 
using the Fourier transform — ^ ^ dkC{k,t)e'^'''^ . We 
get 

-\d^C{k,t) = [k^ + K^] C{k,t) 
Co 

oo 

+ Vn[{k-2nkB)^ + K^]C{k~2nkB,t). (3) 



Considering the weak lattice limit, \vi\ <^ 1, the Bragg 
resonance condition for n = 1 is satisfied for the spatial 
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Fourier modes with peaks centered at fcs and —ks- As- 
suming the resonant initial pulse "^{x, 0) with the Fourier 
indices k G [ks — Afc/2, ks + Afc/2] we get 



: tof [Ci + V1C2] , = Loi [C2 + vlCi] , 

(4) 

where Ci — C{kB+Sk, t) and C2 — C{—kB+Sk, t) are the 
resonant Fourier modes, lji^2 = co[(fcB ± Sk)^ + k^]^ are 
the corresponding frequencies, and 5k € [—Afc/2, Afc/2], 
see Fig. [2] On the other hand, since Eq. (|4]) is the 
second-order system, launching a pulse with the frequen- 
cies lying in the forbidden gap is possible only by sud- 
denly switching on the lattice. The experimental feasi- 
bility of such a setup is an open question, but possible 
in principle: the switch-on time is compared to the Rabi 
period = 27r/r2o with VLq — u}q\vi \ (see Eq. ([6]) below) 
which is much larger than the light period. 
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FIG. 2: (Color online) Bragg resonant coupling of Fourier 
modes (shown by dots). The curve gives the dispersion re- 
lation Lo(k) — cq\/W^+k? (with Co = 0.9 and k = 0.1 in 
dimensionless units). 



Consider first the simplest case with no detuning, i.e. 
the modes with the indices fcs and — fc^. We have 
^{ks) — w(— fcs) = Wo E^nd the system ^ supplemented 
with the initial conditions 

Ci(0) = 1, ^(0) = C2(0) - 0, ^(0) = 0, 

(5) 

corresponding to the initially propagating pulse in a ho- 
mogeneous medium with £ — eq, can be easily solved: 



Ci = cos ( — — 
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where x = ai'g(^'i) — 7''/2 and the terms of order 0(|{>ip) 
are omitted. The Fourier powers of the two peaks (de- 
noted here by |Ci_2p) oscillate with the Rabi frequency 
rio- The oscillations are modulated with the frequency 
2uJo and the amplitude |wi|/4. 

The predictions have been checked by numerical sim- 
ulations of Eq. (m with the truncated expression 
enum(a;) = eo [I — vi cos {2'Kx / d)] (i.e. only the reso- 
nant terms are used). The Gaussian pulse \&(a;,0) = 
explifc^a; — {x — {x))'^/a'^} has been used as the initial 



condition. Fig. [3] shows a good comparison with the 
theoretical solution ([6]) for a pulse with the frequency 
interval Aw « 0.07wo = O.GSrio- 

The sum of the Fourier powers |Cip -I- IC2P = 1 — 
sin(r2ot) sin(2woi) is not conserved since the system 
(ID) with Sk = has the energy 
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|Cip + |C2p+2Re(viCiC2* 



(7) 

For the solution fH]) one can neglect the term 
Re(t)iCiC2) = 0(|uip). Hence, there are oscillations 
between the sum |Cip + IC2P and the analogous sum 
of the time derivatives of Ci^2, with the two frequencies 
LUo i r2o/2, corresponding to the lower and upper band- 
gap ends, and with the amplitude |'0i|/4. The oscilla- 
tions reflect the transfer of the electromagnetic energy 
between the electric and magnetic fields. The Rabi oscil- 
lations can be traced also in the real space, where they 
appear in the form of pulse average position oscillations, 
see Figlllja). The pulse is also spreading (Fig. Hfb)) and 
propagating with a small group velocity. 
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FIG. 3: (Color online) The average Fourier powers Pi, 2 (see 
Eq. (|10|l below) of the numerical solution of Eq. ([l}, the 
thick lines, vs the theoretical solution ([6} with no detuning, 
the thin line. The upper line gives the sum of powers Pi -I-P2. 
Here to = (cofcs)"^, £0 = 1.3, vi = 0.12, and k = Q.05fcs. 
The Gaussian initial pulse is used with the width a — 34d, 
hence AA; = 0.04fcfl, i.e. Alo = 0.65^0. 



As one would expect, for shorter pulses than that of 
Fig. [3] the oscillations show such features as the de- 
phasing and the amplitude damping, since there is a fre- 
quency detuning between the resonant Fourier modes k 
and k — 2kB {k 7^ fcs), see Fig. [51 In the general case, the 
two resonant modes in Eq. ([4]), i.e. Ci{t) — C{kB +Sk,t) 
and C2{t) = C{—kB + Sk, t), have the following form (up 
to the O(|-0i |)-terms) 



Ci = e" 
C2 = e*x-^-7t 



nt\ 7 

cos I — — i — , 
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— , sm — 
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FIG. 4: The pulse average position, defined as 
{x) = J dxx\E\^/ Jdx\E\'^ panel (a), and its half-width, 
[Axf = Jdxix - {x))^\E\^/ Jdx\E\^, panel (b), corre- 
sponding to Fig. |3l 

where 

V 2 +LU^ \vi\ 

(9) 

Since the small 0(|-Ci|)-terms are neglected, the total 
power is conserved in Eq. ^ |Ci p -f IC2 P = 1. Both the 
amplitude of oscillations and the frequency are defined 
by the ratio 7 « ~ of the frequency detuning 

to the band-gap width. Averaging over an interval of 
the size Ak and evaluating the integral by the stationary 
phase method, one obtains for i 3> 1: 

±fcB+Afe/2 

Pi,2 ^ ^ J dA|C(A,t)p 

±feB-Afe/2 

where Auj = uj{kB + Afc/2) - uj{~kB + Afc/2). An 
example of the Rabi oscillations with a large detuning 
(i.e. for a short pulse) is given in Fig. O In this case, the 
average position of the pulse and the pulse width show 
nearly linear dependence on time, the pulse spreads and 
propagates in the lattice significantly. 

In conclusion, an electromagnetic pulse propagation in 
a switched-on photonic lattice with the central frequency 
lying in the forbidden gap is an analog of the Rabi model, 
where the two-level media is the pulse, i.e. its two reso- 
nant Fourier modes, while the photonic lattice serves as 
a classical monochromatic field. The Rabi frequency is 
given by the band-gap width and the oscillations have all 
the characteristic features as in the original Rabi case, 
such as the amplitude damping due to the dephasing. 
An additional feature is the high-frequency modulation 
of the oscillations due to the energy transfer between the 



electric and magnetic fields of the pulse. The effect may 
have applications, for instance, as a base for an all-optical 
trap for light pulses. 
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